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Abstract: Based on a recent formulation of the T^-parameter of a 
photonic crystal fiber we provide numerically based empirical expres- 
sions for this quantity only dependent on the two structural parameters 
— the air hole diameter and the hole-to-hole center spacing. Based on 
the unique relation between the y-parameter and the equivalent mode 
field radius we identify how the parameter space for these fibers is re- 
stricted in order for the fibers to remain single mode while still having 
a guided mode confined to the core region. 
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1 Introduction 

Theoretical descriptions of photonic crystal fibers (PCFs) have traditionally been re- 
stricted to numerical evaluation of Maxwell's equations. In the most general case, a 
plane wave expansion method with periodic boundary conditions is employed [1] while 
other methods, such as the multipole method [2], take advantage of the localized nature 
of the guided modes and to some extend the circular shape of the air-holes. The reason 
for the application of these methods is the relatively complex dielectric cross section of 
a PCF for which rotational symmetry is absent. 

The aim of this work is to provide a set of numerically based empirical expressions 
describing the basic properties such as cutoff and mode-field radius of a PCF based on 
the fundamental geometrical parameters only. 

2 Fiber geometry and numerical method 

Wc consider the fiber structure first studied by Knight et al. [3] and restrict our study 
to fibers that consist of pure silica with a refractive index of 1.444. The air holes of 
diameter d are arranged on a triangular grid with a pitch, A. In the center an air hole 
is omitted creating a central high index defect serving as the fiber core. A schematic 
drawing of such a structure is shown in the inset of the right panel in Fig. 1. 

Depending on the dimensions, the structure comprises both single- and multi-mode 
fibers with large mode area as well as nonlinear fibers. The results presented here cover 
relative air hole sizes, d/A, from 0.2 to 0.9 and normalized wavelengths, A/A, from 
around 0.05 to 2. The modeling is based on the plane-wave expansion method with 
periodic boundary conditions [1]. For the calculations of guided modes presented the 
size of the super cell was 8x8 resolved by 256 x 256 plane waves while for calculations 
on the cladding structure only, the super cell was reduced to a simple cell resolved by 
32 X 32 planes waves. 

3 The y parameter and the relative mode-field radius 

When attempting to establish a simple formalism for the PCF it is natural to strive for a 
result similar to the y-parameter known from standard fibers [4, 5]. However, a simple 
translation is not straight forward since no wavelength-independent core- or cladding 




Fig. 1. left panel shows VpcF calculated from Eq. (1) for d/A ranging from 0.20 
(lowest curve) to 0.70 in steps of 0.05. The dashed line indicates VpcF = i"- The right 
panel shows the relative equivalent mode-field radius, lOpcp/A plotted as function 
of Vpcp for each of the 9 curves in the left panel. The inset shows a schematic 
drawing of the considered PCF structure. 




Fig. 2. The left panel shows curves for constant values of VpcF in a normalized 
wavelength versus relative hole-size plot. The open circles indicate calculated data 
points with full lines to guide the eye. Similarly, the right panel shows curves for 
constant relative equivalent mode-field radius. 

index can be defined. Recently, we instead proposed a formulation of the y-parameter 
for a PCF given by [6] 

VpcF = 27r^y^n|]y(A) -n|gj^(A) (1) 

Although this expression has the same overall mathematical form as known from stan- 
dard fibers, the unique nature of the PCF is taken into account. In Eq. (1), nFM(A) is 
the wavelength dependent effective index of the fundamental mode (FM) and nFSM(A) 
is the corresponding effective index of the first cladding mode in the infinite periodic 
cladding structure often denoted the fundamental space filling mode (FSM) . For a more 
detailed discussion of this expression and its relation to previous work we refer to Ref. [6] 
and references therein. We have recently argued that the higher-order mode cut-off can 
be associated with a value of VpcF = tt [6] and showed that this criterion is indeed iden- 
tical to the single-mode boundary calculated from the multipole method [7]. Recently 
the cut off results have also been confirmed experimentally [8]. Further supporting the 
definition of VpcF is the recent observation [9] that the relative equivalent mode field 
radius of the fundamental mode, wpcf/A as function of T^cf fold over a single curve 
independent of d/A. The mode field radius wpcF is defined as Aqs = mvp^p and cor- 
responds to the width of a Gaussian intensity distribution with the same effective 
area, A^s, as the fundamental mode itself [9]. 

In the left panel of Fig. 1, calculated curves of l^cF as function of A/ A are shown 
for d/A ranging from 0.20 to 0.70 in steps of 0.05. In general, all curves are seen to 
approach constant levels dependent on d/A. The horizontal dashed line indicates the 
single-mode boundary VpcF = tt. In the right panel, wpcf/A is plotted as function of 
VpCF for each of the 9 curves in the left panel and as seen all curves fold over a single 
curve. An empirical expression for wpcF /A can be found in Ref. [9] . The mode is seen 
to expand rapidly for small values of T^pcF and the mode-field radius saturates toward 
a constant value when VpcF becomes large. In fact, it turns out that wpcp/A ~ 1.00 
for VpcF = 1 and wpcf/A ~ 0.62 for Vpcf = tt. In the left panel of Fig. 2, curves 
corresponding to constant values of VpcF are shown in a A/ A versus d/A plot. In the 
right panel, curves of constant wpcp/A is shown, also in a A/A versus d/A plot. Since 
there is a unique relation between ujpcp/A and VpcF [9] the curves naturally have the 
same shape. 

When designing a PCF any combination of d and A is in principle possible. However, 
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Fig. 3. Plot of the parameter space in terms of relative hole size and normalized 
wavelength divided into three regions by the boundaries defined by VpcF = 1 a^nd 
^PCF = In the upper red area the mode penetrates deeply into the cladding 
region and in lower blue region the structure supports a higher-order mode. 



in some cases the guiding will be weak causing the mode to expand beyond the core 
and into the cladding region [10, 11] corresponding to a low value of VpcF- In the other 
extreme, the confinement will be too strong allowing for the guiding of higher-order 
modes [6, 7]. Since both situations are governed by VpcF the design relevant region in 
a A/A versus d/A plot can be defined. This is done in Fig. 3 where the low limit is 
chosen to be VpcF = 1 where wpcp/A ~ 1. How large a mode that can be tolerated 
is of course not unambiguous. However, for lupcF A leakage-loss typically becomes a 
potential problem in PCFs with a finite cladding structure. In non-linear PCFs it is for 
dispersion reasons often advantageous operating the PCF at Vpcp ^ 1 and then a high 
number of air-hole rings is needed to achieve an acceptable level of leakage loss [12]. 

Finally, we note that the practical operational regime is also limited from the low 
wavelength side. In Ref. [13] a low-loss criterion was formulated in terms of the coupling 
length Zc = A/[7ifm(A) — 7ifsm(A)] between the FM and the FSM. In general scattering- 
loss due to longitudinal non-uniformities increases when Zc increases and a PCF with 
a low Zc will in general be more stable compared to one with a larger Zc- Using npM + 
npsM ~ SripM ~ 2nsiiica we can rewrite Eq. (1) as 



from which it is seen that a high value of the IZ-paramcter is preferred over a smaller 
value. In Fig. (3) it is thus preferable to stay close to the single-mode boundary (VpcF ~ 
tt) but in general there is a practical lower limit to the value of A/ A which can be realized 
because when A/A < 0.1 one generally has that Zc ^ A [13]. 

4 parameter expression 

Although the V^-parameter offers a simple way to design a PCF, a limiting factor for 
using Eq. (1) is that a numerical method is still required for obtaining the effective 
indices. In analogy with expressions for standard fibers [5] it would therefore be conve- 
nient to have an alternative expression only dependent on the wavelength. A, and the 
structural parameters d and A. In Fig. 4, we show VpcF as function of A/ A (data are 
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Fig. 4. Plot of VpcF a-s a function of relative wavelength A/A for d/ A ranging from 
0.20 (lowest curve) to 0.80 in steps of 0.05. 

shown by open circles) for d/A ranging from 0.20 to 0.80 in steps of 0.05. Each data set 
in Fig. 4 is fitted to a function of the form 
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and the result is indicated by the full red lines. Eq. (3a) is not based on considerations 
of the physics of the V-parameter but merely obtained by trial and error in order to 
obtain the best representation of calculated data with the lowest possible number of free 
parameters. Prior to the fit, the data sets are truncated at VpcF = 0.5 since wpcF ^ 2A 
in this region (see left panel in Fig. 1) and the data is thus not practically relevant. In 
Eq. (3a) the fitting parameters 3§, and depend on d/A only. In order to extract 
this dependency, suitable functions (again obtained by trial and error) are fitted to the 
data sets for £/, SS, and . We find that the data are well described by the following 
expressions 

w 3.405 x4 , ^ 

^ 0.904 ^ ' 

m ) = 0.200 X f + 0.100 + 0.027 x (1.045 - f )"^'® (3c) 

^ (f)= 0.630 X exp (^^^^I^^ (3d) 

The above set of expressions, Eqs. (3), constitute our empirical expression for the V- 
parameter in a PCF with A/A and d/A being the only input parameters. For A/A < 2 
and VpcF > 0.5 the expression gives values of Vpcp which deviates less than 3% from 
the correct values obtained from Eq. (1). 

5 Endlessly single-mode criterion 

The term endlessly single-mode (ESM) refers to PCFs which regardless of wavelength 
only support the two degenerate polarization states of the fundamental mode [14] . In 
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Fig. 5. Plot of VpcF in the A — > limit as function of the relative air hole size (open 
circles). The full red line represents a fit to the data points and the horizontal 
dash line indicated the ESM limit VpcF = The insert shows a close-up of the 
intersection with the vertical line indicating the air hole size d/A ~ 0.43. 



the framework of the V^-parameter this corresponds to structures for which VpcF < 
for any A/ A [6]. As seen in the left panel of Fig. 1 this corresponds to sufficiently small 
air holes. However, from the plot in Fig. 1 it is quite difficult to determine the exact 
d/A value for which VpcF = tt for A approaching 0. From Eq. (3) it is easily seen that 
the value may be obtained from 



Fig. 5 illustrates this equation graphically where we have extrapolated the data in 
Fig. 4 to A = 0. From the intersection of the full line with the dashed line we find 
that d/A ~ 0.43 bounds the ESM regime. Solving Eq. (4) we get d/k ~ 0.44 and the 
deviation from the numerically obtained value is within the accuracy of the empirical 
expression. 

6 Conclusion 

There are several issues to consider when designing a PCF. In this work we have ad- 
dressed the single/multi-mode issue as well as those related to mode- field radius/field- 
confinement, and mode-spacing. We have shown how these properties can be quantified 
via the 1^-parameter. Based on extensive numerics we have established an empirical 
expression which facilitate an easy evaluation of the T^-parameter with the normalized 
wavelength and hole-size as the only input parameters. We believe that this expression 
provides a major step away from the need of heavy numerical computations in design 
of solid core PCFs with triangular air-hole cladding. 

Acknowledgments 

We thank J. R. Folkenberg for stimulating discussion and M. D. Nielsen acknowledges 
financial support by the Danish Academy of Technical Sciences. 





